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The phenomenon of negative remanent magnetization (NRM) has been observed experimentally in a 
number of heterogeneous magnetic systems and has been considered anomalous. The existence of NRM in 
homogenous magnetic materials is still in debate, mainly due to the lack of compelling support from 
experimental data and a convincing theoretical explanation for its thermodynamic validation. Here we 
resolve the long-existing controversy by presenting experimental evidence and physical justification that 
NRM is real in a prototype homogeneous ferromagnetic nanoparticle, an europium sulfide nanoparticle. We 
provide novel insights into major and minor hysteresis behavior that illuminate the true nature of the 
observed inverted hysteresis and validate its thermodynamic permissibility and, for the first time, present 
counterintuitive magnetic aftereffect behavior that is consistent with the mechanism of magnetization 
reversal, possessing unique capability to identify NRM. The origin and conditions of NRM are explained 
quantitatively via a wasp-waist model, in combination of energy calculations. 



NRM, also known as negative coercivity or inverted hysteresis loop, has been observed in heterogeneous 
thin films, and heterogeneous NP systems (with amorphous matrix, large particle size distribution or 
core shell structure) by various experimental techniques (e.g., Vibrating Sample Magnetometer, SQUID 
magnetometry, Alternating Gradient Force, Magneto-optical Kerr Effect and Hall probe imaging) 1-9 . In normal 
magnetic behavior the magnetization follows the applied field with a phase lag, whereas in magnetic systems that 
display NRM, the magnetization of the descending branch of the hysteresis loop becomes negative while the 
applied field is still positive. Models and theories have been proposed to explain the observed NRM. For example, 
NRM in Ag/Ni multilayer films 4 and granular (Ni, Fe)-Si0 2 films 5 was explained by magnetostatic interaction; 
NRM in YCo 2 /YCo 2 (Fe 10 Ni9o/Fe 10 Ni9o) polycrystal bilayer films and CoNbZr (CoFeMoSiB) amorphous single 
layer films were explained by competing uniaxial (biaxial) and uniaxial anisotropics 10,11 , which was later extended 
to competing cubic and uniaxial anisotropics 1213 ; NRM in Cr 2 0 3 coated Cr0 2 particles 14 , and a randomly 
distributed Co NP system 15 was explained by dipolar interactions between ferromagnetic and superparamagnetic 
particles. 

These previously proposed explanations attributed NRM to heterogeneous structural characteristics, such as 
heterogeneous layers or interfaces, or non-uniform chemical composition of the studied system, and, therefore, 
cannot be applied to a homogeneous nanostructured system. In addition, an agreement on the existence of NRM 
has not been unanimously reached, mainly due to the concern about its perceived violation of the first law of 
thermodynamics, and the lack of adequate rationalization and verification of the inverted hysteresis loop in the 
previous studies. A systematic in-depth investigation into the intriguing NRM phenomenon is still lacking. 

In this paper, we verify the thermodynamic permissibility of the observed NRM, via detailed insight into the 
major hysteresis loops and first-order reversal curves (FORCs) of the studied system and confirm the existence of 
NRM with counterintuitive magnetic aftereffect behavior that is consistent with the magnetic reversal mech- 
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Figure 1 | Structural characterization of 15 nm EuS NPs. (a) HR-TEM micrograph of 15 nm EuS NPs (Insets: histogram of size distribution as fitted into 
function y = — 0.37X 2 + 11.78% — 46.41, and an electron diffraction pattern of face- centered- cubic (fee) EuS NPs, exhibiting features from arising from 
(220) and (200) lattice planes), (b) A high-magnification TEM image of a single EuS NP (Inset: electron diffraction pattern as obtained via Fourier 
analysis) 22 . 



anism. We further propose a novel interpretation of NRM as an 
extreme case of wasp-waist behavior 16 and present a quantitative 
explanation of the origin and conditions of NRM for simple nano- 
magnetic systems with uniform chemical composition and homo- 
geneous structure. A three-dimensionally confined, matrix-free, 
nearly monodispersed EuS NP system is used to verify the presence 
of NRM in a simple homogeneous magnetic system. In contrast to 
direct exchange coupling in 3d ferromagnetic materials, (i.e., nickel, 
iron and cobalt), the magnetic properties of europium monochalco- 
genides (EuX: X = O, S, Se or Te) are caused by the half- filled 4f shell 
of the Eu 2+ ions with a magnetic ground state S 7/2 and are jointly 
governed by ferromagnetic indirect exchange (J x ) among nearest 
neighbor Eu 2+ ions and antiferromagnetic super- exchange (J 2 ) 
among next-nearest neighbor Eu 2+ ions 17,1819 . The behavior of EuX 
nanostructures, including the controversial NRM phenomenon 
observed in EuS NPs, is largely beneficial for the potential magnetic 
and optical applications of EuX crystals on the nanoscale 20 " 22 . 



Results 

Structural characterization. The average size of the EuS NPs, em- 
ployed in this study is 15 nm with less than 15% size distribution. A 
high-resolution transmission electron microscopy (HR-TEM) image 
and a selected area electron diffraction (SAED) pattern of the EuS 
NPs are shown in Fig. 1. 

Magnetic characterization. Zero-field-cooled (ZFC) and field- 
cooled (FC) magnetization measurements from 2 K to 300 K were 
conducted for various field strengths from 20 Oe to 10000 Oe, 
(Fig. 2). When the field was reduced below 1000 Oe, a split be- 
tween the ZFC and FC curves and the formation of a peak, at T p , 
in the ZFC curves were observed. For T < T p , the rotation of spins is 
frozen unless the external applied field is sufficiently large. As the 
temperature is increased in this range, the initial alignment of the 
spins, whose orientations were frozen when cooled to 2 K in the 
absence of an applied field is maintained and the increased 
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Figure 2 | The temperature dependence of the magnetization after zero-field-cooled and field-cooled magnetization using various field strengths from 
20 Oe to 10000 Oe. 
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Figure 3 | Determination of T c . (a) Arrot plots at temperatures between 10 K and 20 K. (b) M s 2 vs. Tplot from which T c for 15 nm EuS NPs was 
determined. T c for bulk EuS is indicated on the graph for comparison. 



thermal fluctuation contributes to the increase in the magnetic 
moment by rotation of some spins that become free to align with 
the applied field. As the temperature increases above T p , all of the 
spins became free to rotate. As a result, the magnetic moments 
decrease as the directions of spins becomes increasingly random- 
ized due to thermal fluctuations. We must emphasize that T p 
cannot be interpreted as a blocking temperature that marks the 
ferromagnetic to superparamagnetic transition, due to the existen- 
ce of hysteresis at temperature significantly above T p 23,24,25 . 

The Curie temperature, T c , of the EuS NPs was determined from 
magnetization isothermals using an Arrott plot 26 , i.e. a M 2 vs. H/M 
plot, as shown in Fig. 3a. The square of the spontaneous magnetiza- 
tion, M s 2 , was determined from the intercept of a linear extrapolation 
from the high field portion (H > 5000 Oe) of the isotherms to H/M 
= 0. A positive (negative) value of the intercept indicates that T<T C 
(T > T c ), i.e., the NPs are in ferromagnetic (paramagnetic) state. For 
the 15 nm EuS NPs, the positive (negative) value of the intercept 
were obtained from the isothermals at T = 16K(T = 17 K); hence 
16K<T C <17K. The obtained M s 2 values were then plotted as a 
function of T (Fig. 3b). The plotted points fall on a straight line with 
coefficient of determination, R 2 = 0.9999, and T c , the temperature at 
which M s 2 = 0, was determined for 15 nm EuS NPs to be 16.4 K, 
which is slightly lower than that of bulk EuS (16.6 K) 17 . Others have 
used isothermal measurements under applied fields up to 5000 Oe to 
evaluate T c of 15 nm EuS NPs 27 , however, our data indicate that 
isotherms measured in the vicinity of T c continue to bend for applied 
fields higher than 5000 Oe, as shown in Fig. 3a. Extrapolations of the 
isothermals measured only up to 5000 Oe would place the evaluated 
T c in an incorrect temperature range, i.e. 15K<T C <16K. Hence, 
when the Arrott plot is applied, one must extrapolate from a suffi- 
ciently high field portion of the isotherms 28 . 

Magnetic hysteresis behavior. Major hysteresis loops were 
measured at temperatures in the range of 2 K to 300 K. The 
temperature dependence of remanent magnetization, M r , is shown 
in Fig. 4. Positive remanent magnetization (PRM) and NRM were 
respectively observed in temperature ranges of 2 K < T < 4.8 K and 
4.8K<T<25K. As the temperature increased from 2 K, the 
magnitude of M r reached zero at the reversal temperature, T r = 
4.8 K, from a positive value, and the temperature at which hystere- 
sis vanished, T hv = 25 K, from a negative value. Zero magnitude of 
M r persisted for temperatures above T hv . Note: T c is in the 
temperature range that exhibits NRM. M r is determined by the 
magnetization of two sublattices that have different temperature 



characteristics. The occurrence of NRM is attributed to sufficient 
antiferromagnetic coupling between the sublattices. 

For T r < T < T hv , the major hysteresis loops were partially 
inverted, with crossover between the descending and the ascending 
branch at a positive field and a negative field. For example, at 10 K, 
crossing between the descending branch and the ascending branch 
can be observed at 30000 Oe and -30000 Oe (Fig. 5). Insets (a) and 
(b) of Fig. 5 respectively show a high field region that displays normal 
hysteresis and a low field region that displays inverted hysteresis. 

FORCs were measured, at various temperatures from 2 K to 
300 K, to reveal elementary processes that cannot be sufficiently 
represented in major hysteresis loops and verify their thermodyn- 
amic permissibility. In FORCs measurements, a class of minor hys- 
teresis loops was obtained by reserving the applied field at various 
points, H r , on the descending (ascending) major hysteresis loop and 
subsequently going to positive (negative) saturation. For T r < T < 
T hv , the temperature range in which NRM is exhibited, the magnet- 
ization along the descending (ascending) branch of the major hys- 
teresis loop is observed to be reversible for different H Y between the 
two crossover points. 

We point out that the major hysteresis loop progressed in a normal 
counterclockwise direction; the occurrence of an apparent clockwise 
major and minor hysteresis loop in the inverted region, that might 
seem to be operating above unity efficiency, i.e. could be used to 
realize perpetual motion machines, violating the first law of ther- 
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Figure 4 | The temperature dependence of the remanent magnetization. 
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Figure 5 | Major hysteresis loop measured at 10 K. Insets: (Left) Detail of a high field region (36500 to 37000 Oe) that displays normal hysteresis. 
(Right) Detail of a low field region ( — 500 to 500 Oe) that displays inverted hysteresis. 



modynamics, is actually due to the intersection of the descending and 
ascending branches of the counterclockwise major hysteresis loop. 
Remarkably, the observed partially inverted hysteresis loop is not real 
hysteresis behavior, but rather the joint behavior of the high-coercive 
surface spins and the low-coercive core spins of EuS NP, which will 
be discussed in detail. The reversibility in minor loops can only be 
observed while the applied field is not sufficiently high to switch the 
high-coercive component of magnetization. Energy over the entire 
loop is conserved; energy gain in inverted hysteresis in the low field 
region is compensated by the energy loss in normal hysteresis in the 
high field region. Hence, the observed partially inverted hysteresis 
loop is thermodynamically allowed, whereas a fully inverted hyster- 
esis loop is not thermodynamically allowed, due to the inevitable 
energy loss in reversing the hard component of magnetization 29 . 

Magnetic aftereffect behavior. We report, for the first time, the 
unique capability of magnetic aftereffect experiments to identify 



the existence of NRM, and present the observation of an unusual 
magnetic aftereffect behavior of EuS NPs in the temperature range of 
T r <T< T hv that is consistent with the mechanism of magnetization 
reversal. A typical magnetic aftereffect behavior of nanomagnetic 
materials that display a slow decay rate 30 is shown in Fig. 6a; the 
intensity of magnetization under a constant applied field decreases 
from its initial value, i.e. the magnetic moment on the descending 
hysteretic curve, towards its anhysteretic ground state with time. 
When the holding field is in the vicinity of the coercive field of the 
descending hysteretic curve, H cd , also known as the negative coercive 
field in typical hysteresis loop, the decay rate is maximized 31,32 and 
the initial magnetization is close to zero. 

If NRM exists, i.e. H cd > 0, at a holding field that is close to H cd , 
the initial magnetization should be smaller in value than its corres- 
ponding anhysteretic ground state, and therefore the magnetic 
moment is expected to increase with time during thermal relaxation, 
as shown in Fig. 7. 




in® ,n « |T<Tr 



Figure 6 | Magnetic aftereffect measurements, (a) Magnetic aftereffect measured for 24 hours at 2 K for various holding fields in the vicinity of H cd . An 
analytical time window (ATW) from —300 s to 24 hr, i.e. ln(t) from 6 to 1 1.4, were used in the measurements. Note: The first —300 s were not included in 
the ATW to guarantee highly stable temperature and field values. In the selected ATW, the decay curves were observed to be quasi-linear and their slopes 
were used to determine the decay coefficient, (b) Comparison of magnetic aftereffect measured at 2 K, 10 K and 30 K for the holding fields of their 
respective H cdy —75 Oe, 14 Oe and 0 Oe, within the selected ATW. 
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Figure 7 | Simulation of the descending and ascending major hysteretic 
curves and anhysteretic curve in the normalized field region of inverted 
hysteresis and the magnetic aftereffect decay curves for the holding fields 
in the vicinity of H c ^. Inset: Detail of the magnetization decay from the 
initial magnetization on the descending hysteretic curve towards the 
anhysteretic curve (indicated by the arrow). 

Figure 6b juxtaposes the magnetic aftereffect of 15 nm EuS NPs 
measured at three representative temperatures, 2 K, 10 K and 30 K, 
in the temperature range of T < T r , T r < T < T hv , and T > T hv , at the 
holding field of their respective H cd , i.e. —75 Oe, 14 Oe and 0 Oe. In 
contrast to the normal magnetic aftereffect behavior observed at 2 K 
(T < T r ), as expected, the magnetic moment increased with time at 
10 K (T r < T < Thv). Therefore, the existence of NRM in homogen- 
eous EuS NP system is confirmed and the consistency between the 
occurrence of NRM and the mechanism of magnetization reversal is 
validated. 

Discussion 

We propose that the NRM is an extreme case of wasp -waist behavior. 
In a typical wasp -waist hysteresis loop, the width of the loop narrows 
as the applied field goes to zero and then opens up 16 . A wasp-waist 
hysteresis loop can be generated if two ferromagnetic materials, one 
of which has a much larger coercivity than the other, are antiferro- 
magnetically coupled. 

Figure 8a shows an extreme case of the wasp-waist behavior. The 
hysteresis loops of the two ferromagnetic materials, with large and 
small coercivity, are respectively denoted as hard loop and soft loop. 
The saturation magnetization obtained in the soft loop, M soft , is 
larger than that obtained in the hard loop, M hard . We assert that 
when the above criteria are met, a wasp -waist hysteresis loop that 
exhibits NRM can be obtained. For given values of M har d and M so ft, 
the shape of the obtained composite hysteresis loop can be adjusted 
to fit experimental data accurately by varying the exchange- coupling 
parameter, which is defined as the dimensionless ratio of the normal- 
ized equivalent field on one material by the normalized magnetic 
moment of the other material 33 . 

The energy per unit volume of the above described magnetic sys- 
tem that exhibits NRM can be given by 

E = K sm 2 p h + K s sin 2 P s0 - M h H cos(P h -6)- M S H(P S - 6) 

(1) 

-JM h M s cos(P h -P s ) 

where the first two terms represent uniaxial anisotropy energy, the 
next two terms represent Zeeman energy, and the last term repre- 
sents the exchange coupling between the hard and soft components. 
K is the anisotropy constant of the hard component, which can be 
determined from 2K/M = H c , and, K s , the anisotropy constant of the 
soft component is assumed to be negligible, i.e. the system has single 



effective uniaxial anisotropy. / is the antiferromagentic exchange 
coupling constant. M h and M s are the magnetization of the hard 
and soft component, fa and fa are the angles between M h and M s 
and the easy axis of the hard component, and fa 0 is the angle between 
M s and the easy axis of the soft component. 0 is the angle between the 
applied field and the easy axis of the hard component. A schematic 
illustration of the wasp-waist system, with the definition of para- 
meters, is shown in Fig. 8b. 

For the EuS NP system, enhanced or modulated physical charac- 
teristics, such as the magnetization 18,27 ' 34 " 41 , Curie/Neel temper- 
ature 18,36 ' 37 ' 39 ' 40 , Verdet constant 41 ' 42 , or absorption/luminescence/ 
birefringence spectra 43 " 46 , are observed compared to their bulk coun- 
terparts, often due to the contribution of uncompensated surface 
spins to the net properties (magnetization, for example) of the NP. 
Indeed, this has been observed when the nanocrystals' average dia- 
meter approaches a critical size 18 ' 35 * 40 - 41 ' 44 ' 49 . As the size of EuS NPs is 
reduced to the nanoscale regime, the surface-to-volume ratio 
increase significantly and a bimodal magnetization state formed by 
surface spins and core spins, under certain temperature and field 
conditions, can be obtained 18 ' 47 " 49 . Such characteristics could lead 
to a variety of nanomaterials applications, which could supplant 
those based on bulk EuS or on other magnetic materials 36 ' 39,41 ' 45 . 
The anisotropy energy of the three-dimensionally confined homo- 
geneous EuS NP system is primarily attributed to the surface atoms, 
i.e., the anisotropy constant of the core atoms is assumed to be 
negligible. At saturation state, the magnetization of core atoms, 
M core , are larger than that of the surface atoms, M sur f ace . In low energy 
regime, i.e. applied field and temperature are low, the core atoms and 
the surface atoms tend to align antiferromagnetically with each 
other. Hence, the EuS NP system meet the criteria of the wasp-waist 
system exhibiting NRM within appropriate window of related para- 
meters, where the surface and core atoms can respectively be viewed 
as the hard and soft magnetization component. 

The magnetization of the EuS NP system, which is the magnet- 
ization component of the net magnetization along the applied field, 
can be given by 

M = M surface cos(P surface -0) + M core cos(fa ore - 0) (2) 

where parameters denoted by subscripts surface and core are equi- 
valent to those denoted by hard and soft in equation (1). 

The magnetization reversal mechanism observed in EuS NPs is 
attributed to bimodal magnetization state formed by surface spins 
and core spins, which is due to evolving balance and changing dom- 
inance among anisotropy energy, Zeeman energy, and aniferromag- 
netic coupling energy under different field and temperature 
conditions 48 . 

When the applied field is sufficiently large, i.e., saturation field, 
both core spins and surface spins align along the applied field 
(Fig. 8c). In this case, Zeeman energy is dominant. While in low 
energy regimes, i.e., low temperature and low field, anisotropy energy 
is dominant and the surface spins can fluctuate along their easy axis. 
However, at very low temperature, i.e. 2K< T < 4.8 K, the anti- 
ferromagnetic coupling between the surface and core atoms are sup- 
pressed, because the rotation of core spins are frozen upon removal 
of the applied field, and mutual alignment of core spins is sustained 
in absence of external applied field. Within appropriate temperature 
window, i.e. 4.8 K < T < 25 K, at the remanent state, core spins align 
antiparallel with the surface spin, i.e. easy axis, and the total magnet- 
ization is negative (Fig. 8d). 

In summary, the NRM observed in a EuS NP system and simple 
nanomagnetic systems with uniform chemical composition and 
homogeneous structure, in general, has been rationalized. A wasp- 
waist model, in combination with energy calculations, is proposed to 
give unambiguous and quantitative explanations to the origin and 
conditions of NRM. 
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Figure 8 | Interpretation of the origin and conditions of NRM. (a) Illustration of an extreme case of wasp-waist behavior that exhibits NRM. 
(b) Schematic illustration of the proposed wasp-waist system. M soft , M hard , and M net are respectively denoted by green, blue and red arrows, (c) Schematic 
illustration of a special case of spin orientation: both surface and core spins align along the field direction. Magnetization of surface and core spins is 
denoted by blue and green arrows, respectively. Total magnetization is denoted by red arrow, (d) Schematic illustration of a special case of spin 
orientation: the core (surface) spins align antiparallel (parallel) with respect to the easy axis and the total magnetization is negative. Magnetization of 
surface and core spins is denoted by blue and green arrows, respectively. Total magnetization is denoted by red arrow. 



Methods 

Sample preparation. Nearly monodispersed EuS NPs were colloidally synthesized 
using a one-step solvothermal synthesis technique 22 ' 50,51 . The synthesis was carried 
out using a Schlenk line. 0.6 mmol of diethylammonium diethyldithiocarbamate, 
0.2 mmol europium oleate, 0.6 mmol of phenanthroline, were dissolved in a solution 
of 1 ml of 1-dodecanethiol and 6 ml of oleylamine contained in a three-neck flask. 
The reaction mixture was then was degassed by purging at 80°C with argon flow for 
45 min. Subsequently, the system was heated to 320°C under vigorous stirring. The 
reaction was allowed to proceed for 1 hour at 320°C, after which the mixture was 
transferred to a glass vial with a glass syringe. 10 mL acetone was added to the fresh 
purple EuS NP sample in oleylamine as synthesized. The sample was precipitated out 
on the bottom of the vial. The upper liquid was subsequently removed, and the EuS 
NPs were then cleaned with acetone for four times. The cleaned EuS NPs were 
dispersed in toluene for electron microscopic measurements. 

Structural characterization. HR-TEM along with SAED pattern measurements were 
conducted using a Philips CM 200 operating at 200 kV. The sample for HR-TEM 
study was made by dropcasting of a EuS NP suspension in toluene onto carbon coated 
copper grids. The two dominant SAED ring patterns respectively correspond to (200) 
lattice plane and (220) lattice plane of FCC EuS NPs. 

Magnetic characterization. Magnetic properties of the EuS NP system were studied 
using a SQUID magnetometer (Quantum Design, MPMS XL 5.0 Tesla). Major 



hysteresis loop and FORCs measurement protocols were adopted in characterization 
of magnetization as a function of applied field. In FORCs measurements, the 
magnetization was measured while the applied field was cycled between various H r on 
the descending (ascending) major hysteresis loop and the positive (negative) 
saturation. ZFC and FC measurement protocols were adopted in characterization of 
magnetization as a function of temperature. In ZFC measurements, the EuS NPs were 
cooled in zero applied field to 2 K, and the desired applied field was applied, i.e. from 
20 Oe to 10000 Oe. Magnetization was measured with ascending temperature to 
300 K. FC measurements were conducted in a similar manner except that the EuS 
NPs were cooled in the presence of an applied field. In magnetic aftereffect 
measurements, the EuS NPs were first saturated by a large positive field, i.e. 50000 Oe, 
which is then rapidly decreased to a field in the vicinity of H cd , where it was held 
constant while the magnetization was measured as a function of time. 

Modeling method. The major hysteresis loops and FORCs that correspond to the 
observed NRM was simulated as an extreme case of wasp-waist behavior using a 
Preisach model. In the Preisach model, the EuS NP was considered to be a 
parallel collection of independent squre-loop hysterons, where each hysteron is 
multiplied by a weight distribution function and then summed 32 . Characteristics 
of the hard and soft components within the Presiach model, i.e. hysterons that 
have high and low coercivities, can be respectively determined upon the 
identification of the Preisach parameters of the measured major hysteresis loops 
and FORCs. 
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